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We review basic constraints on the acceleration of ultra-high-energy (UHE) cosmic 
rays (CRs) in astrophysical sources, namely the geometrical (Hillas) criterion and re- 
strictions from radiation losses in different acceleration regimes. Using the latest avail- 
able astrophysical data, we redraw the Hillas plot and figure out potential UHECR 
accelerators. For the acceleration in central engines of active galactic nuclei, we con- 
strain the maximal UHECR energy for a given black-hole mass. Among active galaxies, 
only the most powerful ones, radio galaxies and blazars, are able to accelerate protons 
to UHE, though acceleration of heavier nuclei is possible in much more abundant lower- 
power Seyfert galaxies. 
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1 Introduction 

The origin of ultra-high-energy (UHE; energy 8 > 10 19 eV) cosmic rays (CRs) remains 
unknown despite decades of intense studies (see e.g. Ref. |1] for a comprehensive re- 
view and Ref. [2] for a recent pedagogical introduction). Recent studies, notably the 
observation of the Greisen-Zatsepin-Kusmin [3JE] cutoff by the HiRes experiment [5], 
further supported by results of the Pierre Auger observatory (PAO) [6j, suggest that at 
least a large fraction of UHECRs is accelerated in cosmologically distant astrophysical 
sources. The birth of (scientific) UHECR astronomy, however, awaits our firm under- 
standing of energies and primary composition of the observed cosmic rays as well as 
identification of at least one astrophysical object where these particles are accelerated. 

Given experimental ambiguities, it is important to understand theoretically, which 
astrophysical objects may serve as UHECR accelerators. It has been understood long 
ago that the UHECR sources should be distinguished objects with extreme physical 
conditions. One simple criterion is the geometrical one: the particle should not leave the 
accelerator before it gains the required energy. Presumably, the particle is accelerated 
by the electric field and confined by the magnetic one; then the geometrical criterion 
is expressed in terms of the particle's Larmor radius which should not exceed the 
linear size of the accelerator (see, e.g., Ref. [7]). In the context of UHECRs, this 
condition is recognized as the Hillas criterion [8] and is often presented graphically 
in terms of the Hillas plot where the accelerator size R and the magnetic field B 
are plotted. We note that even quite recent reviews use either cut-and-pasted or 
slightly refurbished versions of the original 25-years-old plot. However, astrophysics 
experienced enormous progress, if not a revolution, during these decades. One of the 
aims of this study is to give an updated version of the Hillas plot with references to 
either - when possible - measurements or estimates of the magnetic fields and sizes 
of potential astrophysical accelerators. The most important update corresponds to a 
wide variety of active galaxies whose sizes and magnetic fields differ by many orders of 
magnitude from one object to another so that some of them may, while most of them 
may not, accelerate particles to UHE. 

Another restriction on the cosmic-ray accelerators is posed by the radiation losses 
which inevitably accompany the acceleration of a charged particle. The corresponding 
constraints were studied, in particular, in Refs. [U [9j [10] lllj . The radiation losses 
depend on the particular field configuration and the maximal achievable energy of a 
particle in the loss-limited regime depends on the acceleration mechanism. Restricting 
to particular mechanisms or particular field configurations may result in would-be 
contradictory results, cf. Refs. [101 lllj . In this work, we review the radiation- loss 
constraint for different cases; they further limit the acceptable region on the updated 
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Hillas plot. One of possible applications of these general constraints is a study of active 
galaxies correlated with the Auger events [12j . 

One should always keep in mind that even if both geometric and radiation con- 
straints are satisfied, they do not yet guarantee particle acceleration to the correspond- 
ing energy. Each particular source should be discussed in the context of an acceleration 
mechanism operating there. 

The rest of the paper is organised as follows. In Sec. [2j we review constraints on 
potential UHE accelerators, that is model-independent Hillas geometrical constraint 
and limitations due to radiation losses for particular acceleration mechanisms. In 
Sec. [3l we take advantage of the modern astrophysical data and redraw the Hillas 
plot supplemented by the radiation-loss constraints. Our results are summarized and 
discussed in Sec. 2] while brief conclusions are given in Sec. [5j Appendix [A] contains 
derivation of some formulae. 

2 General constraints from geometry and radi- 
ation 

An accelerator of UHECR particles should satisfy several general constraints which 
may be briefly summarized as follows: 

• geometry — the accelerated particle should be kept inside the source while being 
accelerated; 

• power — the source should possess the required amount of energy to give it to 
accelerated particles; 

• radiation losses — the energy lost by a particle for radiation in the accelerating 
field should not exceed the energy gain; 

• interaction losses — the energy lost by a particle in interactions with other 
particles should not exceed the energy gain; 

• emissivity — the total number (density) and power of sources should be able to 
provide the observed UHECR flux; 

• accompanying radiation of photons, neutrinos and low-energy cosmic rays 
should not exceed the observed fluxes, both for a given source and for the diffuse 
background (in particular, the ensemble of sources should reproduce the observed 
cosmic-ray spectrarrQ) . 




The primary concern of this paper is the geometrical and radiation-loss constraints 
(others are briefly quoted when relevant); both restrict the magnetic field and the size 
of the accelerator and can be graphically represented on the Hillas plot. 

2.1 The Hillas criterion 

The Larmor radius Rl of a particle does not exceed the accelerator size, otherwise the 
particle escapes the accelerator and cannot gain energy further. This Hillas criterion 
sets the limit 



for the energy £ gained by a particle with charge q in the region of size R with the 
magnetic field B. 

1 We note that the spectrum of cosmic rays accelerated in a particular source may be very different from 
the spectrum observed at the Earth, cf. Ref. [T5] . 




£ < £n = gBR 



(1) 
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2.2 Radiation losses 



While Eq. ([T]) is a necessary limit, more stringent ones may arise from the energy 
losses: the maximal energy £"i oss a particle can get in an accelerator of infinite size is 
determined by the condition 

rfg (+) = 

dt dt ' 1 J 

where the energy gain rate in the effective electric field E = rjB is (in the particle- 
physics c = 1 units which we use throughout the paper) 

(kept explicit in equations, the efficiency coefficient rj is set to one in figures to obtain 
conservative (optimistic) limits for a given magnetic field - the electric fields in astro- 
physical objects are much less studied observationally compared to the magnetic ones, 
but it is always expected that E <C B). Depending on particular conditions in the 
accelerator, the maximal energy £ max of a particle is limited either by geometrical or 
by energy-loss arguments: 

£max = min{£ H , ^loss} • 

The general expression for total radiation losses for a particle with velocity v moving 
in arbitrary electric E and magnetic B fields reads (see e.g. Ref. [14] ) 

-^-s^M™^ -<*->')■ w 

where q and m are the particle's charge and mass; cross and dot denote vector and 
scalar product, respectively. By making use of relativistic equations of motion, Eq. @ 
can be conveniently rewritten |15j as 



d£(~) 2 - 2 



^F 2 -(F-v) 2 



dt 3 m 2 (l — v 2 

The force F acting on a particle is further decomposed as F = Fii + Fj_, where we 
determine the parallel Fii and perpendicular F^ components with respect to v, that is 
(Fj_ • v) = 0. Then 



dt: • 2 - 



( F l + Fhl-v 2 )). (5) 



dt 3 m 2 (l — v 2 

It is apparent that the second term (the so-called curvature radiation) is suppressed 
with respect to the first one (synchrotron radiation) by an extra power of (1 — v 2 ) and 
therefore may be neglected in the ultrarelativistic regime unless the synchrotron term 
is zero or very small itself. The synchrotron losses are dominant for any generic field 
configuration; however, in a very specific regime v || E || B they vanish, and the losses 
are then determined by the curvature radiation. 



2.3 Different acceleration regimes 

Depending on the scenario of acceleration, we will consider diffusive (stochastic) and 
inductive (one-shot, or direct) mechanisms (see e.g. Ref. [8j for a general discussion of 
these two approaches to UHECR acceleration). 
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The prime examples of diffusive processes are the Fermi first-order [16j and second- 
order (shock, e.g. |17j ) acceleration. Other possibilities include interaction with medium 
by crossing a boundary between layers with different velocities [18] and even transfor- 
mation of a particle into a different one [19| . A recent review and more references can 
be found in Ref. [20]. 

In inductive mechanisms, the particle is accelerated by the large-scale electric field 
continuously and then leaves the accelerator. Strong ordered fields on relatively large 
scales are required; example scenarios are given e.g. in Refs. [211 1221 [231 I2H 122] • For 
our purposes, it is convenient to separate the inductive-acceleration scenarios into two 
groups, depending on whether the configuration of the accelerating field corresponds 
to synchrotron- (e.g. large-scale jets |21| ) or curvature-dominated (neutron stars [22] 
and black holes [231 [21 [25]) losses. 



2.3.1 Diffusive acceleration. 

The losses in this regime are the most serious. This scenario cannot be realized in 
strongly ordered field configurations with v || E || B, therefore the losses are determined 
by the synchrotron limit, 

- ^— = ^-^2 (—) = |-^4^ 2 ^ 2 (synchrotron). (6) 

at 3 R Z L \m ) 3 m 4 

This regime has been studied in Ref. [11] where it has been shown (see Appendix IA.2|) 
that the maximal energy is limited by 

Diffusive mechanisms are quite generic and may work in every realistic environment 
which can host, e.g., a shock wave. Eq. ([7]) does not rely on a particular acceleration 
mechanism and gives a (hardly reachable) upper limit for the maximal energy. 



2.3.2 One-shot acceleration with synchrotron-dominated losses. 

In this regime, the energy loss rate is given by Eq. ([6|) and, given Eq. Q, Eq. ([2]) 
results in the maximal energy 

e - = yi^ 1/2 " 1/2 - <«) 

This acceleration mechanism requires ordered fields throughout the acceleration site; 
its practical realization for UHECR may work in jets of powerful active galaxies [21] . 



2.3.3 One-shot acceleration with curvature-dominated losses. 

The energy loss rate is determined (see Appendix lA.ip by 

d£H 2q 2 (£\ A 

(curvature), (9) 



dt 3 r 2 \ m 

where r is the curvature radius of the field lines which is supposed to be of order of 
the accelerator size and Eq. ([2]) results in the maximal energy 

^=Q) 1/4 ^J^ jR V2 r? V4. (1Q) 
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This mechanism requires ordered fields of very specific configurations which, however, 
may be present in the immediate vicinity of neutron stars and black holes |22[ [2"3[ 1241 



2.4 Summary of results for the maximal energy 

Let us summarize the expressions for the maximal energy £ max (in the comoving frame) 
atteinable by a nuclei with atomic number Z and mass A in the accelerator of size R 
filled with magnetic field B, for different acceleration regimes: 

£h(B,R), B < B (R); 
S loss (B,R), B>B (R), 



£max{B, R) 

where 



B (R) = 3.16 x 10- 3 G (— J fl\ 

is determined from Eqs. dU) and flTJ, © or ([TO]) by requiring £b(B,R) = £i oss (B,R); 
the Hillas constraint is 

S U (B,R) = 9.25 x 10 23 eV Z ( r— ) (77 



kpcy vg 



and the radiation-loss constraints are 

£i oss (B,R) = £ d (B,R) = 2.91 x 10 16 eV'— lf- 



Z 4 Vkpc j \G 
for diffusive acceleration, 

A 2 /B\~ 1/2 , 

£l oss (B,R) = £ S (B,R) = 1.64 x 10 20 eV - ry 1 / 2 

for inductive acceleration with synchrotron-dominated losses and 



1/2 / r>\ 1/4 



f loss (B, 72) = £e(B, i?) = 1.23 x 10 22 eV ^ (J^\ f^j ^ 



for inductive acceleration with curvature-dominated losses. Applications to particular 
objects and graphical representations of the constraints will follow in Sec. U] (see in 
particular Figs. I8l-[j~2]). 

Note that the critical value Bq(R), at which £yl(B, R) = £\ oss (B, R), is the same for 
all three acceleration regimes: in this case, the Larmor radius Rl and the size of the 
accelerator R are equal; within our approximation they coincide also with the curvature 
radius r of the field lines. Therefore the diffusive acceleration regime merges the one- 
shot acceleration because the particle interacts with the shock wave only once in this 
limiting case; moreover, Eqs. © and Q coincide and the two regimes of inductive 
acceleration result in similar losses. 



3 Magnetic fields in particular sources 

A number of astrophysical sources have been proposed where acceleration of cosmic 
rays up to the highest energies can take place (see e.g. Refs. [23 [57] for reviews and 
summary). In this section, we review experimental information on their magnetic 
fields in order to put them in proper places on the Hillas plot. General methods of 
astrophysical magnetic-field studies are discussed e.g. in Refs. [28, 29J; however, a much 
wider variety of them is used for studies of individual sources. 
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3.1 Neutron stars, pulsars and magnetars 



Neutron stars host the highest known magnetic fields in the Universe. In particular, 
magnetars (including anomalous X-ray pulsars) may possess kilometer-scale fields B ~ 
10 14 G and higher while normal neutron stars have B ~ (10 11 -j- 10 12 ) G. Observational 
evidence for these high fields is discussed e.g. in Sec. 6.3 of Ref. [30]. We note also 
a direct (though not widely accepted) observational method to measure B in neutron 
stars, observation of spectral lines giving evidence for resonant Compton scattering at 
cyclotron frequency in the high-field media, see e.g. Ref. [31] for normal neutron stars 
and Ref. [32j for anomalous X-ray pulsars. 

3.2 Active galaxies 

For the purposes of the present study, we will use a simplified classification of active 
galaxies (see textbooks [331 El] and, for a more detailed discussion, Ref. [35] ). Clearly, 
there are many intermediate states and peculiar objects which do not fit this classifica- 
tion well; while they should be studied individually if suspected to be UHECR sources, 
their parameters of relevance (sizes and magnetic fields) are expected to interpolate 
between those of better classified active galaxies. 

Seyfert galaxies - spiral galaxies with bright emission-line nuclei; radio-weak; do 
not possess large-scale relativistic jets; often have starburst activity. 

Radio galaxies - radio- loud elliptical galaxies with relativistic jets. According 
to Ref. J3B], they are classified into two luminosity classes: FR I (less powerful; jets 
brighter towards core; jets may be curved) and FR II (most powerful; straight jets 
brighter at the hot spots at their end points). 

Blazars - (almost) point-like objects with non-thermal spectrum; strongly variable; 
similar in total power to radio galaxies; may be associated with radio galaxies whose 
jets are pointed towards the observer. They may be divided into BL Lac type objects 
(relatively low power; no emission lines; possible counterparts of FR I radio galaxies) 
and optically violently variable quasars (extremely powerful; may have emission lines; 
possible counterparts of FR II) . 

Low-power active galaxies (Seyferts) are much more abundant than radio galaxies 
and blazars. 

Possible acceleration sites in active galaxies include both the central engine (imme- 
diate vicinity of the black hole and the accretion disk) and extended structures (jets, 
lobes, hot spots and jet knots). We will discuss separately the black-hole environment 
and extended structures because of very different conditions for particle acceleration. 
Let us note that the term "active galactic nuclei" (AGN) is often used to describe a 
region much larger than just the black hole and its accretion disc and often includes 
inner jets (or sometimes even larger structures) which we consider separately. 

3.2.1 Supermassive black holes and their environment. 

Measurements of magnetic fields in the central regions of galaxies have been performed 
by means of the following methods (see Fig. Q] for particular results). 

1. Synchrotron self-absorption. Under certain conditions, the low-energy cutoff 
in the spectrum of a compact source may be detected and its shape may be proven to 
correspond to the absorption of synchrotron photons on themselves. If this is the case, 
then the magnetic-field strength may be determined by means of the Slysh formula [13] 
or its modifications. The method works best of all for strong radio sources with resolved 
nuclear components [IT | W2 \ 1231 . 
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Figure 1: The size-field diagram for central regions of active galactic nuclei. Grey colors (box 
and arrows) correspond to Seyfert galaxies, blue colors (error-bar lines marked I) correspond 
to FRI radio galaxies, green colors (error-bar lines marked II) correspond to FRII radio 
galaxies and quasars. Arrows: upper limits from the Zeeman splitting in megamasers (light 
grey, Ref. [37]; medium grey, Ref. [SB]; dark grey, Ref. [39J). Dark blue (I) vertical error 
bar: Faraday rotation measurements, Ref. [40J. The grey box [41J, light blue (I) vertical |42j, 
dark green (II) [43J and diagonal (I and II) jS] error bars correspond to the measurements 
by the synchrotron self-absorbtion method (see text). The allowed region for acceleration of 
10 20 eV protons is located between thick red lines H and S (the lower line, H, corresponds 
to the Hillas limit; the upper one, S, corresponds to the radiation-loss limit for one-shot 
acceleration with synchrotron-dominated losses). 
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2. Polarimetry. Measurements of the Faraday rotation and of resulting depolar- 
ization give estimates of the magnetic field provided the plasma density is known from 
independent observations [3D]. 

3. Zeeman effect in megamasers. Megamasers are compact sources of coherent 
radiation in molecular clouds inside or around the accretion disk. Current precision 
allows to put very stringent constraints on the magnetic fields in these regions from 
non-observation of the Zeeman splitting in megamasers in nearby Seyfert galaxies |37t 
EH EH]. 

4. The iron K a line. Measurements of the width and shape of this X-ray line may 
provide important information about circumnuclear dynamics; in particular, it may be 
used to estimate the magnetic field, though present constraints are quite weak [36] , 

All these direct measurements, however, cannot probe the most interesting region 
in the immediate vicinity of the central black hole, a few Schwarzschield radii (Rs) 
from the center. This region is particularly important because theoretically motivated 
configurations of electric and magnetic fields may allow for negligible synchrotron radi- 
ation of accelerated particles and thus for (relatively weak) curvature-radiation losses. 
Lack of our understanding of the field structure in the accretion disk is transformed 
into uncertainties in the inferred magnetic fields Bbh at the black-hole horizon (see 
e.g. Ref. |47] for a summary of models used for this extrapolation). Direct estimates 
of Bbh are therefore not only scarce but also model-dependent. 

On the other hand, parameters of the environment of a black hole and in particular 
the value of Bbh depend strongly on the black- hole mass Mbh- A conservative upper 
limit on -Bbh follows from the condition that the maximal rate of extraction of the 
rotational energy of a black hole does not exceed the Eddington luminosity [IB] (see 
Ref. [49J for a detailed discussion), 

Bbh < 3.2 x 10 8 (¥$Z\ ^ G. (11) 

Quite old but popular models estimate the Mbh _ Bbh relation from the pressure 
balance (radiation pressure equals to the magnetic- viscosity pressure) JSUJ EI] : 

*"~ 1 *(^f)~ I/ ' G - < 12 > 

An efficient method to constrain the relation between Mbh and Bbh was found in 
Ref. [52J in the frameworks of a particular (not generally accepted) theoretical model 
in which both Mbh and Bbh are related to the observable luminosity at 5100 A. It 
gives somewhat lower values of Bbh than Eq. ()12|) : the best fit is 

log (^p) = (9.26 ± 0.39) - (0.81 ± 0.05) log , (13) 

where the central values of the coefficients are taken from Ref. [52J and the error bars are 
estimated by us from their data. For two cases when rather firm and model-independent 
values of Bbh could be inferred from the observations (synchrotron self absorption 
measured at different radii down to 0.1 pc and extrapolated to Rs, Ref. [44J ) , we 
estimated the corresponding Mbh and found that both values are in a good agreement 
with Eq. (|13p . though precision is quite low. 

Estimates of Bbh versus Mbh are summarized in Fig.0 We will use the upper limit, 
Eq. (jlip . to estimate Bbh for a given Mbh; we note however that realistic values of 
Bbh are 1. . . 2 orders of magnitude lower. Since for the curvature-dominated radiation 
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Figure 2: Magnetic field 5bh at the black-hole horizon versus the black hole mass Mbh- 
Triangles are estimates of Ref. [52] (determined in the frameworks of a particular, not gener- 
ally accepted model) and the grey line (marked [52]) represents their best fit, Eq. ({TBI . Two 
points with error bars correspond to experimental estimates of -Bbh, Ref. [H], using the syn- 
chrotron self-absobtion method (Mbh estimated by us using the stellar velocity dispersion 
from HyperLEDA [53J (thick dark red, FRI radio galaxy 3C 465) and 2MASS K s magnitude 
quoted in NED [M] (thin light green, FRII radio galaxy 3C 111); see Ref. [12] for details). 
Thin blue line (SS) corresponds to the Shakura-Sunyaev estimate, Eq. fll2p . Thick red line 
(Z) represents the Znajek upper limit, Eq. (JTTj) . This conservative upper limit is used in our 
estimates of the maximal cosmic-ray energy. 
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Figure 3: Upper limit on the maximal atteinable energy of protons (red solid line), carbon 
nuclei (blue dotted line) and iron nuclei (orange dashed line) for acceleration with curvature- 
dominated losses near a supermassive black hole, Eq. (I15p . The maximal energy obtained in 
numerical simulations in a particular mechanism [25J is shown by red boxes (protons) and 
orange stars (iron nuclei); these data were obtained from Figs. 5 and 10 of Ref. [25J and 
Eq. (TiTj) of this paper. 



losses higher B always results in higher £ maX ) this assumption is conservative for our 
purposes. 

The size R of the potential acceleration region (that is, the region occupied by 
E || B fields suitable for curvature-dominated losses) is of order R$] therefore both 
R and B are gouverned by a single parameter Mbh> so one may express the maximal 
energy through A/bh using results of Sec. 12.41 Assuming 

R~5R s ~5x 1(T 5 pc A f BH , (14) 

1O 8 M v ; 

one finds that for any reasonable Mbh (ranging from ~ 10 6 M Q for normal galaxies 
through (10 7 . . . 1O 8 )M for Seyfert galaxies to (10 9 . . . 10 10 )Mq for powerful radio 
galaxies and quasars) the maximal energy is determined by radiation losses rather 
than by the Hillas condition and equals to 

£ max = S c c 3.7 x 10 19 eV^u ^ • (15) 



z 1 / 4 v lo8M © 

This general constraint is presented in Fig. for different nuclei (A,Z); for compari- 
son, results of numerical simulations of particle acceleration near a supermassive black 
hole [25] are also plotted. 

In a way similar to other observational manifestations of supermassive black holes, 
both details of cosmic-ray acceleration and radiation losses may depend on the accretion 
rate, accretion mode, environment etc. However, we are interested here in the upper 
limit on the maximal attainable energy of a cosmic-ray particle which is determined 
by Mbh as we have just demonstrated. 
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3.2.2 Jets and outflows of active galaxies. 

Active galactic nuclei fuel large-scale (from sub-parsec to kiloparsec and even Mega- 
parsec length) extended more or less linear jets. Revolutionary progress in the angular 
resolution of radio (sub-milliarcsecond) and X-ray (sub-arcsecond) imaging allowed for 
detailed studies and modelling of physical conditions in jets. We briefly review basics 
of the current understanding of jet properties following Refs. [55, 56J. 

Seyfert galaxies possess extended structures which are often non-collimated (open- 
ing angle of 45° or more) and are found to be non-relativistic; they are sometimes 
determined as "outflows" , reserving the term "jets" to strongly collimated relativistic 
flows. X-ray emission from these outflows, when present, is well described by thermal 
radiation (sometimes associated with star- forming regions in the outflow [57J). 

Relativistic jets reveal themselves in non-thermal X-ray emission studied now in 
great detail. The jets are spatially resolved into components; in nearby jets (Cen A) 
inner and outer layers and bright knots are resolved. It is often assumed that all jets 
are fuelled by the central black hole; the energy flux is dominated by the magnetic- 
field energy at sub-parsec scales but becomes particle-dominated at parsec scales. The 
emission of low-luminosity sources (FR I radio galaxies and BL Lacs) is adequately 
described by the synchrotron models from radio to X rays; their jets are decelerated by 
entrainment of gas and dissipate in the end. High-power FR II and quasar jets bring 
their energy flux directly to their terminal hot spots and require additional (e.g. Comp- 
ton) component to describe their spectra. Comparison of radio to X-ray observations 
gives rather firm evidence to the origin of the emission of FR I jets from accelerated 
particles and to acceleration of these particles not only in a finite number of shocks 
but also by means of some distributed mechanism along the jet [581 159] . Quite rarely, 
relativistic jets are present in exceptionally powerful Seyfert galaxies; in these cases, 
they have properties very similar to FR I jets [56]. Models of multifrequency spectra 
allow to constrain the magnetic field, the key parameter of the synchrotron radiation. 
The estimates depend also on the electron density; the degeneracy is often removed 
either by the equipartition assumption or by a simultaneous measurement of the self- 
Compton component, when applicable. When error bars are given, they include the 
corresponding uncertainties. Some of these estimates [5SI HEB EH [62j [63] are presented 
in Fig. SI In some cases, existence of ordered fields through the jet was proven, so that 
the inductive acceleration may be possible (sec e.g. Rcf. |21j). 

3.2.3 Jet knots, hot spots and lobes of powerful active galaxies. 

When a relativistic jet is present, it may be accompanied by internal shock regions 
(knots), terminal shock regions (hot spots) and extended regions in the intergalactic 
space fuelled by the jet after its termination (lobes). These regions are typically absent 
in low-power active galaxies (Seyfert galaxies): knots are observed mostly in jets of 
FR I radio galaxies and quasars, lobes are typical for radio galaxies, hot spots are 
present in the most powerful FR II radio galaxies and quasars. Magnetic fields may be 
determined either by X-ray synchrotron observations alone (assuming equipartition) 
or by combined multifrequency observations of both synchrotron and Compton radia- 
tions (allowing to relax the equipartition assumption which occurs, in the end, a good 
approximation, see e.g. Ref. [64j ffi A summary of measurements [66 } \67\ [68] is given 
in Fig. [51 

2 An interesting approach to determination of the magnetic field in a knot in M87 65 exploits the energy 
dependence of the energy loss rate, assuming it is synchrotron-dominated. The resulting ~ 0.6 mG field is 
in a good agreement with equipartition-based estimates. 
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Figure 4: The size-field diagram for jets and outflows of individual active galaxies. Grey 
colors correspond to Seyfert galaxies, data from Refs. [56] (light-grey diagonal line), [60] 
(grey vertical error bar), [6T] (short dark grey diagonal). Blue box corresponds to FR I radio 
galaxy [62]; green triangles represent quasar jets [63]. The allowed region for acceleration of 
10 20 eV protons is located between thick red lines H and S (the lower line, H, corresponds 
to the Hillas limit; the upper one, S, corresponds to the radiation-loss limit for inductive 
acceleration with synchrotron-dominated losses). 
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Figure 5: The size-field diagram for knots (triangles), hot spots (boxes) and lobes (circles) of 
individual powerful active galaxies. Filled black symbols correspond to quasars and blazars, 
filled grey symbols correspond to FR II radio galaxies, empty symbols correspond to FR I 
radio galaxies (data from Ref. [66J, X-ray observations assuming equipartition); boxes with 
error bars represent the "best-guess" estimates of Ref. [HZ]; dotted grey circles correspond 
to FR II lobes studied in Ref. [BS] (comparison of radio and X-ray observations without the 
equipartition assumption). The allowed region for acceleration of 10 20 eV protons is located 
between thick red lines H and D (the lower line, H, corresponds to the Hillas limit; the upper 
one, D, corresponds to the radiation-loss limit for diffusive acceleration). 
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Figure 6: The size-field diagram for Galactic star-forming regions (R < pc) and starburst 
galaxies (R > 0.1 kpc). Thick (error-bar) lines correspond to measurements of the Zeeman 
splitting in masers (grey, Ref. [69J; red, Ref. [70]; blue, Ref. [TTj; green, Ref. [72\). Black 
dots represent results of submillimeter imaging polarimetry of Ref. [75]. Data for normal 
(empty boxes), starburst (triangles) and extreme starburst (filled boxes) galaxies are taken 
from Ref. [75]; black symbols correspond to the minimal-energy field estimates while grey 
symbols correspond to equipartition field estimates. The allowed region for acceleration of 
10 20 eV protons is located between thick red lines H and D (the lower line, H, corresponds 
to the Hillas limit; the upper one, D, corresponds to the radiation-loss limit for diffusive 
acceleration). 

3.3 Star formation regions and starburst galaxies 

Measurements of the magnetic field in Galactic star-forming regions becomes possible 
with the Zeeman splitting in masers in circumstellar disks \69\ [TP] [7T] 172] and infrared 
imaging polarimetry [73]. Though these regions in our Galaxy have never been con- 
sidered as possible sites of UHECR acceleration, these measurements may give some 
hints to the fields in larger star-forming regions in starburst galaxies, where particles 
could be accelerated to very high energies e.g. in shocks from subsequent supernova 
explosions [73]; magnetic fields in these extragalactic sites are measured indirectly. A 
summary of measurements is given in Fig. [6j a number of arguments in favour of higher 
(equipartition) fields in starburst galaxies were presented in Ref. [75] while continu- 
ity (see Fig. [6]) with the Galactic measurements may support lower (minimal-energy) 
estimates. 

3.4 Gamma-ray bursts 

Estimates of the magnetic field in gamma-ray bursts (GRBs) assume [76] that the ori- 
gin of both prompt and afterglow emissions in a certain part of the spectrum is the 
synchrotron radiation of relativistic electrons. This assumption is supported by mea- 
surements of the afterglow spectra and lightcurves and by observation of the strongly 
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Figure 7: The size-field diagram for gamma-ray bursts. Horizonthal lines represent estimates 
of Ref. [76] which assume the synchrotron origin for the prompt emission (dark grey) and 
the afterglow (light grey). The allowed region for acceleration of 10 20 eV protons is located 
between (the lower lines, H, correspond to the Hillas limit; the upper ones, D, correspond to 
the radiation- loss limit for diffusive acceleration). Dashed lines assume T = 500, thick lines 
assume T = 50, thin lines assume T = 1. 

polarized prompt emission (see Ref. [76] for discussion and references). Ref. [76] 
quotes B ~ 10 6 G for R ~ (10 13 -j- 10 15 ) cm (prompt emission) and B ~ 1 G for 
R ~ (10 16 4- 10 18 ) cm (afterglow) (we assume that the estimates correspond to the 
observer's rest frame). Another, somewhat higher field estimate may be obtained fol- 
lowing Ref. [10] (see also Ref. |77| ) from the total luminosity of a GRB, assuming that 
the magnetic- field energy 8 m is a fraction e m < 1 of the radiation energy £ TSu d- However, 
this estimate depends strongly on the assumed beaming. 

Within the scope of this paper, we may estimate the maximal energy f max of ac- 
celerated particles in the comoving frame following equations of Sec. 12.41 for shock 
(diffusive) acceleration. The GRB shells are however ultrarelativistic (r ~ 100, see e.g. 
Ref. |78j) and we have to multiply the comoving- frame £ max by T to get the maximal 
rest-frame energy. Results are presented in Fig. [7J which, for the GRB case, is more 
instructive than the summary plots of Sec.[H We note that at large T, the maximal en- 
ergy may be limited by interactions with thermal photon field (not taken into account 
in the present work) and decreases as at large F |10j . 

3.5 Galaxy clusters, superclusters and voids 

Information about cluster magnetic fields comes mostly from observations of their ex- 
tended radio, and sometimes X-ray, emission. These observations are reviewed e.g. in 
Refs. [281 1791 180] ■ where more references to original works may be found. Estimates 
based on equipartition (see e.g. Refs. [Hl[82]), as we ^ as those assuming Compton scat- 
tering on CMB photons, favour values of B ~ (0.1 -j- 1) /J.G at megaparsec scales; Fara- 
day rotation measurements (see e.g. Refs. [831 i84l [85]) favour somewhat higher fields, 



16 



B ~ (1 -T- 5) fiG. Model-dependent numerical simulations remain the main source of 
information about magnetic fields at the supercluster scales (R ~ 100 Mpc), escepially 
in voids. Estimates vary between B ~ 10" 11 G [86] and B ~ 10~ 8 G [87]. 

4 Summary and discussion 

Based on the data collected in Sec. [3] and on the limits on the maximal energy, Sec. 12. 4[ 
we redraw here the Hillas plot supplemented by radiation-loss constraints. Figures [8]- 
[TOl give constraints for particular acceleration regimes while Figs. [11] [12] represent our 
updated summary Hillas plots. 

The weakest possible constraints (for inductive acceleration with curvature-dominated 
losses) are presented in Fig. [8] Constraints for inductive acceleration with synchrotron- 
dominated losses, applicable mostly to inner and outer jets of active galaxies, are given 
in Fig. El while constraints for the most general diffusive acceleration are presented in 
Fig. Ell Figure [11] represents our version of the Hillas plot with constraints for 10 20 eV 
protons. Figure [T2l is the same plot but for 10 20 eV iron nuclei. 

Constraints for neutron stars follow from Sec. 13. 1\ even for the least restrictive 
acceleration regime they are not satisfied for UHE particles. In active galaxies, various 
regimes of acceleration may operate. In the immediate vicinity of the central black 
hole (up to a few Rs), the field configuration allows for the inductive acceleration with 
curvature-dominated losses. These regions are denoted as "BH" in Figs. [8] [TTT fT2| the 
parameters correspond to Eqns. Q3] and [11] The latter one is an upper limit on the 
field, so we extend the boxes for two orders of magnitude lower in B, cf. Fig. [2] Beyond 
a few Rs, the E || B field structure no longer holds, but coherent fields may still be 
present in inner jets. For these central parsecs of AGN (denoted as "AD" in Figs. F9 H12 j) 
we use the field estimates from Fig. [1] For the extended parts of active galaxies (jets, 
jet knots, hot spots and lobes) we use field estimates summarized in Figs. H] The 
summary boxes for starburst galaxies include both equipartition and minimal-energy 
estimates (Fig. [6]). For GRB, the summary plots present synchrotron-based estimates 
for both inner and outer shocks (see Fig. [JJ for a more instructive plot). Field estimates 
for clusters, superclusters and voids follow Sec. 13.51 

The constraints discussed here and expressed in terms of the Hillas plot are neces- 
sary, but they should be supplemented by other limits (listed in the beginning of Sec. [2]). 
We note that in estimation of the maximal atteinable energy, important constraints 
are put by interactions of accelerated particles with ambient photons. In particular, 
interaction with the cosmic microwave background is important for large, R >Mpc, 
sources (lobes of radio galaxies, clusters and voids), while interaction with the internal 
source radiation field is important for ultraluminuous sources (GRB and AGN). These 
constraints, considered elsewhere, further restrict the number of potential UHECR ac- 
celerators^. For the diffusive shock acceleration, these constraints have been studied 
e.g. in Ref. [9]. 

The maximal energy for the supermassive black holes is readily expressed in terms of 
a single parameter, the black-hole mass Mbh- We used the upper limit on the magnetic 
field, i?BH; which is most likely one or two orders of magnitude higher than the actual 
values, so that our estimate, Eq. (fT5|) . is robust. It depends weakly (y R/Rg) on the 
assumed size of the acceleration region. 

While we tried to make all constraints as robust as possible, it is clear that they 

3 In certain cases the proton-gamma interactions, instead of pure dissipation, can significantly amplify the 
acceleration process [E5] ■ 
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Figure 8: The size-field plot with constraints from geometry and radiation losses for the 
regime where losses are dominated by curvature radiation. These are minimal possible 
losses and these constraints are therefore the most weak. Boxes denote parameter regions 
for objects in which conditions for this loss regime may be satisfied, that is immediate 
neighbourhood of neutron stars (NS), anomalous X-ray pulsars and magnetars (AXP) and 
of supermassive central black holes (BH) of active galactic nuclei, from low-power Seyfert 
galaxies (Sy) to powerful radio galaxies (RG) and blazars (BL). The shaded area corresponds 
to the parameter region where acceleration of protons to 10 20 eV is possible. Lines bind from 
below the allowed regions for 10 19 eV protons (thin full line), 10 20 eV iron nuclei (thick full 
line), 10 18 eV protons (dashed line) and 10 17 eV protons (dotted line). Right-hand parts 
of the lines represent the Hillas constraint while left-hand (steeper) parts represent the 
radiation-loss constraint. 
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Synchrotron-dominated losses 
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Figure 9: The size-field plot with constraints from geometry and radiation losses for the 
regime of one-shot acceleration with synchrotron-dominated losses. Boxes denote parameter 
regions for objects in which conditions for this loss regime may be satisfied, that is central 
parsecs (AD) of active galaxies (low-power Seyfert galaxies (Sy) and powerful radio galaxies 
(RG) and blazars (BL)) and relativistic jets of powerful active galaxies. The shaded area 
corresponds to the parameter region where acceleration of protons to 10 20 eV is possible. 
The line binds the allowed regions for 10 20 eV iron nuclei. Lower lines represent the Hillas 
constraint while upper (horizonthal) lines represent the radiation- loss constraint. All quan- 
tities are given in the comoving frame, so the maximal energy for jets should be multiplied 
by the bulk Lorentz factor of the jet which may be as large as ~ 10 for leptonic jets and 
~ 100 for hadronic jets [TO] . 
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Figure 10: The size-field plot with constraints from geometry and radiation losses for the 
regime of diffusive acceleration with synchrotron-dominated losses. Boxes denote parameter 
regions for objects in which conditions for this loss regime may be satisfied, that is central 
parsecs (AD) of active galaxies (low-power Seyfert galaxies (Sy) and powerful radio galaxies 
(RG) and blazars (BL)), relativistic jets, knots (K), hot spots (HS) and lobes (L) of power- 
ful active galaxies (RG and BL); non- relativistic jets of low-power galaxies (Sy); starburst 
galaxies; gamma-ray bursts (GRB); galaxy clusters and intercluster voids. The shaded area 
corresponds to the parameter region where acceleration of protons to 10 20 eV is possible. 
The line binds the allowed regions for 10 20 eV iron nuclei. Lower lines represent the Hillas 
constraint while upper lines represent the radiation- loss constraint. All quantities are given 
in the comoving frame, so the maximal energy for jets and shells of GRBs should be multi- 
plied by the bulk Lorentz factor which may be as large as ~ 10 for leptonic jets and ~ 100 
for hadronic jets and GRBs [TP] . 
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1(T U eV protons 
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Figure 11: The Hillas plot with constraints from geometry and radiation losses for 10 20 eV 
protons. The thick line represents the lower boundary of the area allowed by the Hillas 
criterion. Shaded areas are allowed by the radiation-loss constraints as well: light grey 
corresponds to one-shot acceleration in curvature-dominated regime only; grey allows also 
for one-shot acceleration in synchrotron-dominated regime; dark grey allows for both one- 
shot and diffusive (e.g. shock) acceleration. See captions to Figs. El [10] for notation of 
boxes corresponding to potential sources. 
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Figure 12: The same as Fig. [TTJbut for 10 20 eV iron nuclei. The most important difference 
with Fig. [TT] is that acceleration of iron nuclei to 10 20 eV is possible (unlike for protons) in 
low-power active galaxies (e.g. Seyfert galaxies). 
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should be considered as order-of- magnitude estimates (in fact, typical precision of the 
magnetic-field determination is an order of magnitude) and, for individual unusual field 
configurations, can be quantitatively violated. An example of such a configuration is a 
linear accelerator with curvature radius r of field lines exceeding the size of the source 
R; then R should be substituted by r in Eq. ([9|). The estimates should be used with 
care for the cases when the magnetic field changes violently within the accelerator (for 
instance when particles are accelerated by magnetic reconnection). A more detailed 
modelling of acceleration and losses is required in these cases. 

One of our most important conclusions is that low-power active galaxies (e.g. Seyfert 
galaxies) cannot accelerate protons to energies > 5 x 10 19 eV. Indeed, when the ex- 
tended structures (jets and outflows) are present, the magnetic field there is far too 
weak to satisfy the Hillas condition (even for very rare relativistic jets), see Fig. HI The 
same is true for the accretion disks, where the field is nicely constrained from above by 
non-observation of the Zeeman splitting in megamasers (Fig. [1]) . The most favourable 
conditions for acceleration correspond to the immediate vicinity (a few Rs) of the 
central black hole where the upper limit on the maximal energy is given by Eq. (|15p . 
Since for Seyfert galaxies, Mbh i$ (10 -r- 1O 8 )M , proton acceleration to ~ 10 20 eV is 
not allowed. However, these (and only these) central parts of Seyfert galaxies can, in 
principle, accelerate protons to ~ 10 18 eV and heavy nuclei to ~ 10 20 eV, if interactions 
with ambient photons are weak enough. Though heavy nuclei are much less abundant 
than protons, Seyfert galaxies themselves are much more abundant, and hence typically 
close to the observer, than powerful radio galaxies and blazars, so that their population 
can contribute to the UHECR spectrum. 

5 Conclusions 

We reviewed constraints on astrophysical UHE accelerators and presented the Hillas 
plot supplemented with radiation-loss constraints and updated with recent astrophysi- 
cal data. Contrary to previous studies, we emphasised that active galaxies span a large 
area on the plot, and only the most powerful ones (radio galaxies, quasars and BL Lac 
type objects) are capable of acceleration of protons to UHE. If UHECR particles are 
accelerated close to the supermassive black holes in AGN, then most likely the mech- 
anism is "one-shot" with energy losses dominated by the curvature radiation. Other 
potential UHE acceleration sites are jets, lobes, knots and hot spots of powerful active 
galaxies, starburst galaxies and shocks in galaxy clusters. Acceleration of particles 
in supercluster-scale shocks, gamma-ray bursts and inner part of AGN is subject to 
additional constraints from pj interactions which are not discussed here. 

Unlike protons, heavy nuclei can be accelerated to UHE in circumnuclear regions 
of low-power active galaxies. Since these galaxies are abundant, this contribution to 
UHECR flux may be important, leading to a mixed primary cosmic-ray composition 
at highest energies. 

The authors are indebted to D. Gorbunov, S. Gureev, A.M. Hillas, V. Lukash, 
A. Neronov, S. Popov and D. Semikoz for interesting discussions. We acknowledge the 
use of online tools \53\ [54"] . This work was supported in part by the grants RFBR 07- 
02-00820 09-07-08388 (ST), NS-1616.2008.2 (ST) and by FASI government contracts 
02.740.11.0244 (ST) and 02.740.11.5092 (KP, ST) and by the Dynasty Foundation (KP, 
ST). . 
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A Derivation of electrodynamic results 

A.l Energy losses for the curvature radiation (see e.g. 
Ref. [BE!) 

Consider a particle moving along curved field lines. The particle has a longitudinal 
velocity component (vii || B) and a drift component (vj _L B). This drift component 
provokes the appearance of the Lorentz force which curves the particle's trajectory 
towards the field lines. For a relativistic particle, 



v\\ 2 m ( £ 



Vd 



qBr \m 



so the Lorentz force is 



F L = q[v d x B], 

v 2 m / f 
r \m / 

The energy losses are in general determined by Eq. ([5]) which may be rewritten as 

— = F — (F ■ v) 

dt 3m 2 (l-v 2 ) I 1 1 

In the regime we consider, (F • v) = and consequently 

dt 3r 2 \ m 
In the ultrarelativistic limit v\\ — > c one obtains Eq. 



A. 2 The maximal energy for diffusive acceleration [11 



Consider a flow propagating through a magnetized medium. An accelerated particle 
gains energy by repeated scattering off the flow. After every scattering, the particle 
travels along the Larmor orbit, radiates and slows down according to Eq.([6]); conse- 
quently 



f d£ 2," I R 2 



hence 

1 = 1 + ^1 [ R B 2 ( x) d X 
The maximal energy £ = £ CI is determined by setting £q — > oo, 

and we obtain Eq. ([7J). 
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